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Abstract 

In this paper, a new geometry index of Very Long Baseline Interferometry (VLBI) 
observing networks, the volume of network V, is examined as an indicator of the 
errors in the Earth orientation parameters (EOF) obtained from VLBI observations. 
It has been shown that both EOF precision and accuracy can be well described by the 
power law a = aV^ in a wide range of the network size from domestic to global VLBI 
networks. In other words, as the network volume grows, the EOF errors become smaller 
following a power law. This should be taken into account for a proper comparison of 
EOF estimates obtained from different VLBI networks. Thus performing correct EOF 
comparison allows us to accurately investigate finer factors affecting the EOF errors. 
In particular, it was found that the dependence of the EOF precision and accuracy 
on the recording data rate can also be described by a power law. One important 
conclusion is that the EOF accuracy depends primarily on the network geometry and 
to lesser extent on other factors, such as recording mode and data rate and scheduling 
parameters, whereas these factors have stronger impact on the EOF precision. 

1 Introduction 

Very Long Baseline Interferometry (VLBI) is a unique technique to determine the Earth ori- 
entation parameters (EOF), because it is the only method, which can provide the complete 
set of five EOF: terrestrial pole coordinates Xp and Yp, Universal Time UTl, and celestial 
pole coordinates Xc and Yc, with highest accuracy. Therefore, investigation of various influ- 
ences on the errors in EOF derived from the VLBI observations is important for improving 
the results. One of the main factors affecting the EOF results is the network geometry. 

Evaluation of the dependence of the EOF errors on the network geometry is important 
for many practical tasks, such as comparison of EOF estimates obtained on different VLBI 
networks, observation scheduling, design of VLBI experiments and networks, optimal use of 
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the existing VLBI stations for determination of EOP. For this reason, an analysis of the dif- 
ferences in EOP values and errors obtained from different observing programs and networks 
is one of the most actual topics in geodetic VLBI, and many authors have contributed to 
these studies, e.g. MacMillan and Ma (2000); Sokolskaya and Skurikhina (2000); Johnson 
(2004); Lambert and Gontier (2006); Searle (2006). However the authors of these papers 
did not investigate directly the immediate connection between EOP and the quantitative 
parameters of the network geometry. 

As follows from the basic principles of VLBI measurements, errors in geodetic and astro- 
nomical parameters derived from the VLBI observations directly depend on the number and 
the length of the baselines involved. Let us have a look at the basic VLBI equation, keeping 
the main terms in the geometric delay: 

cr = f ■ (Qel) = r (e; ■ (Qe;)) , (1) 

where r = re^ is the baseline vector in terrestrial coordinate system, r and el- are its length, 
and corresponding unit vector, is a unit geocentric vector in celestial coordinate system 
in the direction of the radio source, Q is the rotation matrix for the transformation from 
celestial to terrestrial coordinate system. To get a corresponding equation of condition to 
be used in a parameter solution, one has to compute the partial derivatives of Eq. ([T]) with 
respect to the parameters to be solved. In case of EOP, the partial derivatives are given as 

where p is one of the Earth rotation parameters. 

As one can see from Eq. ([2]), the partial derivatives of the interferometric delay with 
respect to EOP are proportional to the baseline length. This means that a longer baseline 
yields smaller EOP errors. It should be mentioned, however, that since we consider the VLBI 
observations carried out on the real near-spherical Earth, the longest baselines compatible 
with the Earth's diameter in practice do not necessarily have smallest EOP error due to the 
problem of the common visibility of radio sources from the two ends of the baseline. For this 
reason, an optimal baseline length can be found for every specific geodetic or astronomical 
task, but this analysis is out of the scope of this paper. 

To perform a more detailed analysis of the dependence of the EOP errors on the baseline 
orientation, one can rewrite Eq. ([1]) as 

CT = Te cos 6 cosh + rpSinS , (3) 

where and are the equatorial and polar components (projections) of the baseline vec- 
tor r respectively, 6, h being the source declination and hour angle respectively. Detailed 
consideration of this equation which can be found e.g. in Dermanis and Mueller (1978), 
Ma (1978), Nothnagel et al. (1994), Schuh (2000), along with a joint analysis of the set of 
Eqs. ([3]) written for all baselines of the given station network, allows us to conclude that 
larger baseline equatorial component yields smaller UTl error, and larger both polar and 
equatorial projections are needed to get a better separation and less errors in Xp, Yp, Xc, 
and Yc. Also important is that errors in the pole coordinates Xp and Yp may depend on the 
network longitudinal orientation in case of a moderate longitude span of the network. 
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In general, one can expect that the sought-for dependence of the EOP errors on the 
network geometry is a function of several indices, such as: 

• the network span in various directions e.g. geodetic (latitude A(y9, AA) or Cartesian 
(AX, AY, AZ) coordinates; 

• the network orientation, e.g. central longitude, prevailing baselines directions; 

• the number of stations or baselines. 

All these factors definitely influence, to a greater or lesser extent, the EOP results obtained 
from VLBI observations. Hence they should be accounted for in order to provide more 
rigorous comparison of the EOP derived from different networks and separate the impact 
of the network geometry and other factors, such as recording data rate and scheduling 
parameters, on the EOP errors. However, a large number of the network parameters makes 
it difficult to develop a practical method of comparison of EOP. Therefore, it would be useful 
and convenient to have a generalized index of the network geometry which could serve as an 
argument of a simple enough but sufficiently accurate function describing the dependence of 
the EOP errors on the network geometry, and could be used in the comparative studies. 

In this paper, continuing the previous study (Malkin 2007), such a generalized index of 
the VLBI network geometry is considered, namely the volume of network, more precisely 
the volume of the polyhedron with the network stations at the vertices. In section |2l the 
dependence of the EOP errors on the network volume has been investigated, and in section [3l 
this dependence will be applied to a more rigorous comparison of the EOP errors obtained 
from different VLBI networks. 

It is important to distinguish between two kind of errors — precision, which is a measure 
of the repeatability (reproducibility) of estimates, and accuracy, which refers to the agree- 
ment between our estimates and the true value (veracity). Hereafter, it is considered that 
the EOP precision can be characterized by the formal error (uncertainty) of the EOP esti- 
mates computed making use of a data processing software. To assess the EOP accuracy one 
can compare the EOP obtained from VLBI observations with an independent EOP series 
provided by another space geodesy technique. In this paper, like in other similar studies 
(e.g. Vennebusch et al. 2007), we used for comparison the EOP series provided by the 
International GNSS Service (IGS, Dow et al. 2005). The term 'error' is used here as a 
collective name for precision and accuracy in the case when both precision and accuracy are 
considered. 

For completeness, mention may be made of another measure of the EOP quality, the 
Allan deviation, which may be used either in its original formulation (Gambis 2002) or with 
extensions developed in Malkin (2008) for analyzing unequally weighted and multidimen- 
sional observations, e.g. both pole coordinates simultaneously. This method was not used 
in this work, although it may be worth trying for supplement analysis. 

For our computations, the VLBI observations were used, collected on the global and 
regional networks in the framework of observing programs coordinated by the International 
VLBI Service for Geodesy and Astrometry (IVS, Schliiter and Behrend 2007), and stored in 
the IVS Data Center^]. 

^http://ivscc. gsfc.nasa.gov/products-data/data. html 
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Table 1: The volume of the smallest and largest IVS networks 
Session Volume, Mm^ 



JADE-0610 


9.925E-04 


EUROPE-35 


5.285E-03 


JADE-0601 


1.008E-02 


EUROPE-56 


1.761E-02 


T2037 


4.607E+02 


T2038 


4.723E+02 


T2041 


4.816E+02 


T2043 


4.871E+02 



2 Dependence of the EOP errors on the network vol- 
ume 

This study begins with the computation of an about 10.7-year EOP series for the period 
from July 1996 till February 2007, 1440 24-hour sessions in total. The start epoch of the 
investigated time series was chosen to be the same as the start epoch of the IGS EOP series 
igs00p03.erp used hereafter for accuracy assessment of the VLBI terrestrial pole coordinates. 

During the session processing the volume of the session network was computed simulta- 
neously with EOP in the following way. 

1. Compute the tetrahedron mesh for the network polyhedron by means of the Delaunay 
triangulation making use of the GEOMPACK package by B. Joe (Joe, 1991). 

2. Compute the volume of each tetrahedron as a scalar triple product: 

\{r2 - n) ■ {{rs - fl) X (rl - fl))\ / 6 , (4) 
where rl, 1^,1^3, fl are the geocentric station vectors. 

3. Compute the total network volume as the sum of the volumes of all the tetrahedrons. 

Some examples of the volume of different IVS observing networks processed in this work 
are presented in Tabled where the smallest and the largest IVS networks are shown. The first 
column of the table contains the IVS session name. The smallest IVS networks are typical 
of Japanese (JADE) and European (EUROPE) regional sessions. The largest networks were 
observed in T2 experiments primarily aimed at the improvement of the terrestrial reference 
system. One can see from Table [1] that the volume of IVS networks differ from each other 
by several orders of magnitude. For comparison, the volume of the Earth is 1083 Mm^. 

Now, if we plot the EOP precision as a function of the network volume, as shown in 
Fig. [H, a clear linear (in log-log scale) dependence between these values can be seen. 

To investigate the dependence of the EOP errors on the network volume in detail the 
following computations were performed. 

1. The EOP series mentioned above was arranged in 9 groups by the network volume V 
in Mm^: 1/<0.1, {VTO)''<V<{VTO)''+\ k = -2,...,5; 
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Figure 1: Dependence of the EOF precision on the network volume V. Downward, average 
formal error in the terrestrial pole coordinates Xp and Yp, Universal Time UTl, and celestial 
pole coordinates Xc and Yc are depicted (logarithmic scales on both axes) 
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2. For the sessions fallen into each group the following quantities were calculated: 

• the average network volume; 

• the average formal errors in Xp, Yp, UTl, and Yc regarded as EOF precision; 

• weighted root-mean-square (WRMS) of the differences in Xp, Yp between the 
computed EOF and the IGS combined series igsOOpOS.erpH after removing the 
constant bias between the VLBI and IGS EOF series regarded as EOF accuracy. 

Details and results of computation are shown in Table [2] and Figs. [2] and [31 The nine 
averaged points thus obtained were then used to compute the parameters of a power law 
fitting best to the data 

(T = aV% (5) 

where a is an EOF error, i.e. can stand for precision or accuracy. For actual computation 
this dependence was transformed to the linear form 

logcr = 6 + c logl^ , (6) 

where h = log a. The parameters b and c were computed by means of the least square linear 
fit with weights dependent on number of sessions fallen into each group. The results are 
presented in Tables [3] and HI and depicted as the regression lines in Figs. |2] and [3l One can 
see that the dependence of both EOF precision and accuracy on the network volume can be 
nicely described by a power law. 

Generally speaking, it is natural that the EOF precision improves with larger network 
size (volume) since a large network has larger projections of the network baselines on the 
rotation axes. More interesting and less evident is that the EOF accuracy follows the same 
power law with practically the same exponent c but with different factor a. 

3 Comparison of observing programs 

In this section, an application of found dependence of the EOF errors on network geometry 
to a comparison of EOF obtained from different IVS observing programs will be considered. 

To investigate the dependence of the EOF errors on the observing program, we have 
performed the same computations as described in the previous section, with the only dif- 
ference that EOF results were arranged in 13 groups by the IVS observing programs: 11 
global networks Rl, R4, RD, RDV, NEOS-A, CORE-A, CORE-B, CONT02, CONT05, T2, 
E3, and 2 regional networks EURO, JADE. Although the two last programs are primarily 
intended for geodesy researches in relatively small geographic regions, and are not supposed 
to obtain scientifically useful EOF results, they are included in this study for the following 
reasons: 

• including the relatively small regional networks in this study provides significantly 
larger range of the network volume, and thus allows us to get more reliable estimates 
for the slope of the log-log regression line corresponding to the power law discussed 
above; 

ft p : / /cddis . gsfc . nasa. gov / gps /products/ 
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Table 2: Results of splitting of the VLBI EOP series into 9 groups by the network volume V. 

For each V range, number of sessions N is given, and for each parameter, minimal, maximal 



. average vahies 


are shown dow: 


nward. 


Unit for X^. 


Yp. X,. 


. Yf.: mas 


. imit for UTl: 


V range, Mm^ 


N 


V, Mm^ 




EOP precision 




EOP 


accuracy 








Xp 




UTl 


Xc 


Yc 


Xp 


Yp 


V< 0.1 


20 


0.001 
0.079 
0.037 


0.505 
5.213 
1.599 


0.661 
3.230 
1.500 


0.335 
2.231 
1.109 


0.131 
0.629 
0.383 


0.128 
0.780 
0.404 


3.118 


3.033 


0.1<\/<0.316 


7 


0.134 
0.276 
0.209 


0.115 
0.889 
0.516 


0.068 
1.415 
0.667 


0.020 
0.564 
0.268 


0.056 
0.346 
0.209 


0.055 
0.315 
0.190 


0.415 


0.385 


0.316<y<l 


21 


0.333 
0.992 
0.667 


0.138 
3.451 
0.608 


0.096 
3.735 
0.686 


0.029 
1.495 
0.286 


0.080 
0.910 
0.208 


0.080 
0.979 
0.207 


0.726 


0.491 


1<V<3.16 


13 


1.01 
2.40 
1.48 


0.149 
1.229 
0.431 


0.119 
1.599 
0.482 


0.051 
0.691 
0.189 


0.084 
0.617 
0.191 


0.085 
0.487 
0.180 


0.481 


0.482 


3.16<y<10 


59 


3.33 
9.24 

6.98 


0.052 
1.380 

0.232 


0.067 
1.282 

0.228 


0.018 
1.363 

0.195 


0.031 
0.267 

0.095 


0.025 
0.272 

0.094 


0.251 


0.237 


10<V<31.6 


188 


10.6 
31.2 
22.1 


0.052 
0.835 
0.136 


0.047 
0.627 
0.120 


0.013 
0.744 
0.059 


0.030 
0.303 
0.075 


0.032 
0.353 
0.076 


0.221 


0.188 


31.6<y<100 


698 


32.0 

99.9 
65.4 


0.026 

2.355 
0.114 


0.026 

0.905 
0.094 


0.010 

0.890 
0.049 


0.017 

1.000 
0.065 


0.017 

0.995 
0.065 


0.186 


0.173 


100<V^<316 


422 


101 
313 
170 


0.021 
0.678 
0.072 


0.020 
0.402 
0.069 


0.009 
0.642 
0.035 


0.015 
0.285 
0.049 


0.015 
0.287 
0.049 


0.116 


0.129 


316<\/<1000 


12 


325 
487 
399 


0.030 
0.091 
0.055 


0.034 
0.110 
0.061 


0.015 
0.058 
0.037 


0.027 
0.079 
0.047 


0.027 
0.079 
0.043 


0.106 


0.080 



Table 3: Rep resentation of the EOP precision by a power law log a — b + c log V 



EOP 


Unit 


b 


c 


Xp 


mas 


-0.340 ± 0.022 


-0.351 ±0.018 


Y, 


mas 


-0.342 ± 0.023 


-0.373 ± 0.019 


UTl 


ms 


-0.0603 ±0.0043 


-0.0382 ±0.0036 




mas 


-0.771 ±0.016 


-0.238 ±0.013 


Yc 


mas 


-0.772 ±0.016 


-0.238 ±0.013 
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Figure 2: Dependence of the EOF precision (mean formal errors for the 9 groups of volume) 
on the network volume V (logarithmic scales on both axes). Solid line corresponds to the 
power law (Table [3]) g 



0.01 



0.01 



0.1 



10 



100 



1000 



V, MnT" 




0.01 I — — — — — I 

0.01 0.1 1 10 100 1000 

V, Mm^ 

Figure 3: Dependence of the EOF accuracy (WRMS of the differences with respect to IGS 
EOF for the 9 groups of volume) on the network volume V (logarithmic scales on both axes). 
Solid line corresponds to the power law (Table H]) 



Table 4: Representation of the EOF accuracy by a power law log a = b + c log V 



EOF 


Unit 


b 


c 


Xp 


mas 


-0.212 ±0.045 


-0.315 ±0.038 


Yp 


mas 


-0.263 ± 0.049 


-0.298 ±0.041 
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Table 5: Characteristics of IVS observing programs. Only the programs and observations 
used in this paper are shown. In some columns, the average value is given in parentheses 



Program 


Observation 


Number 


Number 


Number of 


Network volume, Data rate. 




period 


of sessions 


of stations 


observations 


Mm^ 


Mbps 


Rl 


2002-2007 


261 


4-8 (6.2) 


700-4800 (2000) 


2.4-320 (140) 


256 


R4 


2002-2007 


258 


4-8 (6.3) 


400-3200 (1300) 


7.7-273 (84) 


56-128 


RD 


2004-2006 


18 


5-7 (6.1) 


1200-3100 (2300) 


16-27 (22) 


1024 


RDV 


1997-2007 


55 


14-20 (17.5) 7600-29300 (17500) 


28-313 (160) 


128 


NEOS-A 


1996-2001 


283 


4-6 (5.1) 


300-2000 (1100) 


4.7-117 (54) 


56 


CORE-A 


1997-2000 


79 


4-6 (5.5) 


400-2000 (1300) 


1.0-187 (94) 


56 


CORE-B 


1997-2001 


53 


4-8 (6.2) 


400-3500 (1600) 


0.3-235 (76) 


56 


CONT02 


2002 


15 


7-8 (7.9) 


2300-3500 (3000) 


91-92 (92) 


256 


CONT05 


2005 


15 


10-11 (10.9) 


5600-6500 (6000) 


263-265 (265) 


256 


T2 


2002-2006 


41 


5-15 (9.0) 


300-5600 (1800) 


11-487 (219) 


64-128 


E3 


2002-2006 


31 


4-6 (4.5) 


200-800 (400) 


1.0-88 (25) 


128 


EURO 


1996-2006 


36 


4-9 (6.6) 


400-5400 (2800) 


0.01-1.4 (0.6) 


64-128 


JADE 


1999-2006 


12 


4-6 (5.1) 


700-2100 (1400) 


0.001-0.05 (0.03) 


128 



• some IVS analysis centers include EURO, and sometimes JD sessions in their EOP 
series. 

The results shown in the previous section show that taking the small networks into consid- 
eration does not corrupt the estimates of the parameters of the power law describing the 
dependence of the EOP errors on the network volume, which can be seen in Figs. |2] and [31 
A detailed description of different IVS observing programs can be found at the IVS 
websit^, and their main relevant characteristics are shown in Table \5\ Note that in fur- 
ther computation the average value of the extremes was taken for observing programs with 
changing data rate. 

For all the observing programs listed above, except RD (Research & Development ex- 
periments), all the sessions observed during the used time span were included in this study. 
As to the RD program, only 18 sessions observed with the recording data rate 1 Gbps were 
used. It is worth mentioning that CONT02 and CONT05 programs were special 2-week 
continuous observing campaigns primarily aimed at the highest EOP quality. It can also be 
mentioned that all observing programs were observed with the Mark IV terminal developed 
in the USA (Whitney 1998), except E3 sessions observed with the S2 terminal developed in 
Canada (Wietfeldt 1996), and JADE sessions observed with the K4 terminal developed in 
Japan and compatible with Mark IV (Kiuchi 1997). An overview and comparison of VLBI 
terminals can be found e.g. in Petrachenko (2000). The S2 terminal provides maximum 
recording rate of 128 Mbps, whereas the other terminals use the recording rate of 256 Mbps 
for regular experiments and up to 1 Gbps in RD sessions, as can be seen in Table [51 

The results of computations are presented in Figs. [H and [5l where the solid lines corre- 

^http://ivscc. gsfc.nasa.gov/program/descrip. html 
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spond to the power law with the parameters given in Tables [3] and H] for EOP precision and 
accuracy, respectively. Comparing these results with those presented in the previous section, 
one can see that the scatter of the points with respect to the regression line in Figs. H] and 
is greater than that in Figs. [2] and [3], which can be explained by the fact that, as a rule, 
networks of different size participated in the same observing program, which can be clearly 
seen from Table [51 Besides, the difference in some scheduling options, such as optimization, 
recording mode and source selection also influence the EOP errors and cause their deviation 
from the common (average) law. 

Collating the results obtained above for the EOP precision and accuracy one can see that 
the scatter (WRMS of the residuals) of the EOP accuracy with respect to the regression 
line (0.033 mas) is much less than that of the EOP precision (0.054 mas). The difference 
becomes even larger if the WRMS is scaled to the magnitude of the accuracy and precision, 
respectively, more precisely to the a factor in Eq. (JSj) — 0.039 vs. 0.108. This may mean 
that the EOP accuracy follows more strongly a power law than the EOP precision. In other 
words, one can conclude that the EOP accuracy depends primarily on the network geometry 
and to less extent on other factors, which have stronger impact on the EOP precision. 

An analysis of the data presented in Table [5] and Figs. H] and [5] allows us to make some 
interesting perceptions. Here are two examples. 

• The known fact of smaller errors of the EOP obtained from the Rl observations with 
respect to the R4 program can be mainly explained by the difference between the 
average volume of the Rl and R4 networks, 140 and 84 Mm^, respectively. 

• The observations made in the framework of the E3 program, using S2 terminal, in- 
volving relatively small number of stations, and delivering relatively small number of 
observations, show relatively bad precision. However, the accuracy of the E3 EOP, 
after accounting for the network size, is at a level of the accuracy of other observ- 
ing programs involving more stations, using higher recording data rate, and collecting 
much more observations. 

It is interesting to see how the general law derived from the whole data set works for sep- 
arate IVS observing programs. For this comparison, three observing programs were selected, 
Rl, R4, and NEOS-A, each having a large number of sessions and a large range of network 
volume. All the three programs are designed for the highly accurate EOP determination 
and use Mark IV terminals, with different recording data rate however, as can be seen from 
Table [51 Results of computations are presented in Fig. [61 

In these plots, one can see that the three compared observing programs show similar 
slope in the dependence of the formal errors on the network volume, and this slope is close 
to the general power law derived from the whole data set. The small discrepancies in the 
slope between the separate observing programs among them and with the general slope given 
in Table [31 can be evidently explained by the difference in the number of observations and 
the range of the network size. Indeed, finer dependencies of the EOP errors on the network 
geometry may contribute to these discrepancies. One can also see in Fig. [6] that session 
points for R4 and NEOS-A are laying practically on the general regression line, whereas Rl 
sessions are shifted down with respect to the general law, which corresponds to the averaged 
program points in Fig. [H One of the possible reasons of a better precision of the Rl EOP 



11 



R1 

R4 

RD 

RDV 

NEOS-A 

CORE-A 

CORE-B 

CONT02 

CONT05 

E3 

EURO 
JADE 
T2 



0.01 0.1 



10 



100 



V, Mm-' 



1000 



10 



0.1 r- 



R1 


□ 


R4 


■ 


RD 




RDV 


• 


NEOS-A 


A 


CORE-A 




CORE-B 




CONT02 


▼ 


CONT05 




E3 


♦ 


EURO 




JADE 


• 


T2 





0.01 



0.01 0.1 



0.1 



0.001 



0.01 0.1 



1 10 

V, Mm^ 



100 



R1 

R4 

RD 

RDV 

NEOS-A 

CORE-A 

CORE-B 

CONT02 

CONT05 

E3 

EURO 
JADE 
T2 



1000 



R1 


□ 


R4 


■ 


RD 


O 


RDV 


• 


NEOS-A 




CORE-A 


A 


CORE-B 


V 


CONT02 


T 


CONT05 


o 


E3 


♦ 


EURO 




JADE 


• 


T2 





0.01 



0.01 0.1 



0.1 



R1 


□ 


R4 


■ 


RD 


O 


RDV 


• 


NEOS-A 


A 


CORE-A 


A 


CORE-B 




CONT02 


T 


CONT05 


O 


E3 


♦ 


EURO 




JADE 


• 


T2 





0.01 0.1 



10 



V, Mm-' 



100 



1000 



Figure 4: Dependence of the EOF precision on the average network volume V for different 
observing programs (logarithmic scales on bO)2h axes). Solid line corresponds to the power 
law (Table [3]) 
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Figure 5: Dependence of the EOP accuracy on the average network volume V for different 
observing programs (logarithmic scales on both axes). Solid line corresponds to the power 
law (Table iD 
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Figure 6: Dependence of the EOF precision on the average network volume V for three 
IVS observing programs Rl, R4 and NEOS-A (logarithmic scales on both axes). Solid line 
corresponds to the power law (Table [3]) 
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Table 6: Original (cr) and modified (am) errors in pole coordinates (average of Xp and Yp). 

Unit: mas 

Program Precision, Accuracy, 

mas mas 







cr 






cr 




Rl 


0. 


,061 


0, 


.446 


0. 


115 


0, 


.770 


R4 


0. 


,088 


0, 


.520 


0. 


155 


0, 


.853 


RD 


0. 


,079 


0, 


.272 


0. 


170 


0, 


.564 


RDV 


0. 


,046 


0, 


.352 


0. 


137 


0, 


.961 


NEOS-A 


0. 


,105 


0, 


.512 


0. 


178 


0, 


.816 


CORE-A 


0, 


,099 


0, 


.609 


0. 


179 


1, 


.027 


CORE-B 


0. 


,121 


0, 


.691 


0. 


,188 


0, 


.998 


CONT02 


0. 


,051 


0, 


.311 


0. 


111 


0, 


.630 


CONT05 


0. 


,026 


0, 


.243 


0. 


079 


0, 


.672 


T2 


0. 


,118 


1, 


.021 


0. 


,148 


1, 


.178 


E3 


0. 


,253 


0, 


.929 


0. 


223 


0, 


.775 


EURO 


0. 


,691 


0, 


.535 


0. 


,885 


0, 


.693 


JADE 


1, 


,771 


0, 


.435 


4. 


,884 


1, 


.267 



may be a higher recording data rate of 256 Mbps as compared with 56-128 Mbps used during 
NEOS-A and R4 sessions. However, possible dependence of the EOP errors on the recording 
data rate is definitely a finer effect than the dependence on the network volume. 

Thus, it is clear from the results presented above that to correctly compare the errors 
in EOP obtained from different VLBI networks it is necessary to properly account for the 
network size. This can be done by the following method. Before comparison of different 
observing programs, corresponding EOP errors should be reduced to the unit volume. To 
do so, the errors are to be modified as follows: 

am = a/V% (7) 

where a is the original error (precision or accuracy), and a-m is the modified error. Original 
and modified errors in pole coordinates (average for Xp and Yp) computed using Eq. ([7]) 
are shown in Table [61 These results give a quantitative confirmation of the preliminary 
conclusions made from Figs. H] and O 

Comparison of the data presented in Tables and E] allows us to make a conclusion on 
the dependence of the EOP errors on the recording data rate, which is depicted in Fig. [3 
From this data, one can see that the dependence of the errors in pole coordinates on the 
recording data rate can also be represented by the power law: log cr = (0.437 ± 0.041) — 
(0.352 — ±0.121) logR for the pole coordinates precision and logo" = (0.338 ± 0.023) — 
(0.194 — ±0.068) logi? for the accuracy, where R is the recording data rate in Mbps. Again, 
like the case of dependence of the EOP errors on the network volume, the EOP accuracy 
shows stronger dependence on the registration data rate than the EOP precision. This can 
be seen from comparison of the ratio between the largest and smallest modified values of 
precision and accuracy. 
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Figure 7: Dependence of the modified EOP precision (top) and accuracy (bottom) on the 
recording data rate (logarithmic scales on both axes) 
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4 Conclusion 



The volume of a VLBI network, i.e. the volume of a polyhedron formed by the network 
stations, can serve as a generalized index of network geometry. A comparison of EOP 
obtained on different observing networks allowed us to show that both EOP precision (formal 
errors) and accuracy (WRMS with IGS EOP series) strongly follow a power law with respect 
to the network volume in a wide range of the network size, from regional to global one. Thus, 
the network volume seems to be a governing factor of the errors in the VLBI EOP. 

Indeed, the dependence found in this paper cannot predict errors in EOP derived from 
a particular VLBI session, which depend on many reasons, such as details of the network 
geometry and baseline orientation, scheduling parameters (optimization, recording mode and 
data rate, source selection), station operation quality, etc., and may differ from the computed 
value by a factor of 2-3, as can be seen from Fig. [H It is intended for overall comparison of 
the observing networks. 

Based on this result, the method is proposed to correctly compare the EOP errors, both 
precision and accuracy, for different VLBI networks and observing programs. Following this 
method, the errors should be compared after a reduction to the unit volume in accordance 
with Eq. (I7j). After that, comparison of the modified EOP errors obtained from observations 
on different networks (observing programs) does not depend on the network size anymore, 
which allows us to estimate the influence of other factors on EOP results. In particular, 
dependence of the EOP errors on the recording data rate was derived from an analysis of 
the observations made in the framework of several IVS observing programs. It was found 
that this dependence can also be described by a power law. 

It is remarkable that dependence of the EOP accuracy on both the network volume and 
the recording data rate is described by a power law more accurately than that for the EOP 
precision. This may show evidence that the EOP accuracy depends mainly on the network 
geometry, and the EOP precision depends to greater extent on other factors mentioned 
above. 

The results obtained in this paper can also be used for planning geodetic VLBI cam- 
paigns, optimal use of the existing IVS network and for designing new VLBI networks, such 
as VLBI2010 (Petrachenko et al. 2004) and GGOS (Global Geodetic Observing System, 
Pearlman et al. 2006, 2007). 

Of course, the proposed method of accounting for the network geometry is applicable 
only for the networks comprising four and more stations, but this does not seem to be a 
serious limitation since the most valuable scientific VLBI results are obtained on ramified 
networks. 
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